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We investigated the microscopic and electronic structures amorphous oxide semiconductors 
^ InGaZn04 (a-IGZO) and the role of O-deficiency through the first-principle calculations. The 

structure of the amorphous oxide is complicated by the admixture of many different kinds of sub- 

(N 

structures, however it is surprisingly found that the band tail states, which are well-known to 
be present in the amorphous semiconductors, are few generated for the conduction band mini- 
mum (CBM). The electronic structure around CBM is little affected by the disorder and also by 
the O-deficiency. Free electron carriers can be generated without a creation of donor-level in the 
03 O-deficient amorphous oxide. 
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Recently the amorphous semiconductor metal-oxides (ASMO) such as InGaZnO (a- 
IGZO), ZnSnO, and InSnO(ITO) [1-3] are extensively investigated, since they can be base 
materials for the thin-film transistor (TFT) of flat panel display as well as transparent 
conducting oxides (TCO). The n-type conductivity of ASMOs can be much higher than 
that of the conventional semiconductors such clS clS ct - Si :H, since the conduction path of 
electron carriers is characterized mainly by non-directional spherical cation- s orbitals of the 
metal elements [4] . The O-deficiency is suggested to be a important source of high electrical 
conductivity of ASMO and TCO. The carrier concentration in ASMO can be controlled over 
a wide range by the O-deficiency [5]. However the microscopic understanding of the role of 
the O-deficiency is not yet resolved with many unresolved aspects. The optical measure- 
ments indicate that the O-deficiency induces deep levels in ASMO [2], and the first-principles 
calculations also showed that the O-vacancies in a-IGZO can create various deep or shallow 
levels [6]. It is suggested that the deep trap centers by the O-deficiency adversely affects the 
bias and illumination stress instability of oxide TFT [3, 7]. A puzzling aspect is that the 
carrier mobility is measured to be much higher in the disordered a-InGaZnO (a-IGZO) than 
in the well-ordered crystalline IGZO (c-IGZO), unless the O-deficiency is too serious [4]. 
This is not understood by a simple defect picture of O-vacancy, since the cation-driven 
dangling bonds around O-vacancies can generate localized levels which can play a role of 
trap center. Since the microscopic understanding for the O-deficiency is limited, the results 
from the studies of the O-vacancy in ZnO or other crystalline oxides have been applied to 
understand the role of O-deficiency. Therefore, the systematic study about the role of the 
O-deficiency in the ASMO is needed. 

In this Letter, we report on the results of first-principles calculations carried out to 
understand the role of the O-deficiency in a-IGZO as a representative of the ASMO. Our 
approach is different from a previous study in that we investigate the properties of an O- 
deficient a-IGZO structures that are generated via the optimization through a molecular 
dynamics simulated annealing processes. We find that (i) the disordered a-IGZO is mixed 
by many different kinds of building blocks, however (ii) the band tail states, which is well- 
known to be present in amorphous materials, can be surprisingly absent for CBM in the 
amorphous-IGZO, and (iii) the localization character of the CBM is little affected by the 
O-deficiency in a-oxide, and no localized donor-center is formed, that is fundamentally in 
contrast to the O-vacancy in crystalline IGZO (c-IGZO), and (iii) the result of the disorder 
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or O-deficiency is to enhance the formation of localized states around the valence band. We 
propose that these results provide a new perspective to understand the role of O-deficiency 
in ASMO. 

The first-principles calculations were performed by the projector augmented wave (PAW) 
method [8] of the Vienna ab-initio simulation package (VASP) [9]. The Perdew-Burke- 
Ernzerhof exchange-correlation functional (PBE) [10] approach utilizing the generalized 
gradient approximation (GGA) scheme was employed and the LDA+U method was used to 
describe the localized semi-core states of Zn-3d orbitals [11]. These semi-core orbitals are 
not accurately described by normal LDA calculation, in which their overlap with the O-p or- 
bitals is enhanced due to a strong self-interaction effect in the localized orbitals. The strong 
Coulomb energy at the localized d orbitals can be compensated by the LDA+U method. 
We used U=5 eV for only the Zn-3d orbitals. As an accuracy test, the lattice constants 
of the binary oxides, In 2 03, Ga20"3, and ZnO, are calculated to be 12.42 A, 10.117 A, and 
3.283 A, respectively, in a good agreement with the experimental values of 12.23 A, 10.117 
A, and 3.25 A. 

In order to generate a reliable structure for the stoichiometric and O-deficient a-IGZO, 
we repeated the simulated annealing (SA) steps using molecular dynamic calculations on 
supercells. The simulated temperature was increased up to 4000 K, and decreased slowly to 
K and the SA was repeated between 1000 K and K until the total energy and the atomic 
structure were converged. We tested several supercells containing 56, 84, and 112 atoms. 
Calculated results on 112 atoms supercell are given here, since the cell size dependence is 
weak between 84 and 112 atoms cells. The total energy of the a-IGZO generated in this 
manner is higher by only 0.18 eV/atom than that of c-IGZO. 

In order to understand the microscopic structures of a-IGZO, we examined the short-range 
orders by calculating the effective coordination numbers (CN*) of cations and the averaged 
effective bond-length (d*)[12, 13]. We find that the structure of a-IGZO is complicated by 
the admixture of many kinds of C-Oqtss molecular building blocks, where C denotes a cation 
and CN is the coordination number of a cation, as shown in Fig. la. The nondirectional 
cation-s-obitals make these various substructures accommodated in the ASMO. The average 
of d* increases as we go from Ga to Zn to In. In addition, d* increases with increasing CN*, 
and this is due to the increased repulsive interaction between O atoms. The graph of d* 
vs. CN* can be nicely fit by a straight line. Fig. lb show that the average (<CN* >) 
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and the standard deviation of CN* (AstdCN*) increases, as we go from Zn to Ga and to 
In. The In atoms are noted to incorporate the most diverse set of substructures, while Zn 
prefers mainly a tetrahedral geometry. These indicate that the substructure of In is flexible. 
Although In and Ga atoms are isoelectric, the In-<CN* > is larger than Ga-<CN* >, due 
to the Ins larger atomic size. The column-Ill cations donating three electrons can attract 
more O 2- ions than the column-II Zn which donates two electrons. The overall (<CN* >) 
are smaller in a-IGZO than in the c-IGZO. Thus the volume of a-IGZO is larger by 0.6 %. 

Before we examine the O-deficient state of a-IGZO, we investigated the properties of 
an metastable O-vacancy structures, that are formed by the removal of an O atom from 
the stoichiometric structure of a-IGZO. The atomic positions were relaxed through the 
minimization of Hellman-Feynman forces. This case has been tested previously by Kamiya 
et al. [6]. The formation enthalpies of the generated O- vacancies are estimated to be between 
0.77 eV and 4.61 eV(assuming contact with a reservoir of O2 gas). The total energy are 
affected by the cations around V(O). Generally as the number of In atoms surrounding 
the V(O) is larger, the formation energy is lowered, and the more stable vacancies tend to 
have the shallower donor levels, although generally the O- vacancies make donor-like sub-gap 
levels, as suggested by Kamiya[6]. 

Here, we note that the O-deficient a-IGZO is not simply described by the O-vacancy 
structure. These vacancies are metastable states, thus we tested the additional relaxation 
of a-IGZO with a missing oxygen atom in a supercell by repeating the SA process. The 
total energy becomes smaller by 0.45 eV than the most stable V(O) structure. Below, we 
will discriminate between the metastable V(O) structure and the latter O-deficiency state. 
The low formation enthalpy (0.32 eV) of the neutral shallow-donor-like O-deficiency state 
indicates that the (2+)-charged state can be more easily formed, when Fermi level is around 
mid-gap in this wide-gap materials. 

We now examine the structural change by the O-deficiency in a-IGZO, with respect to 
the short-range order. We find, as shown in Fig. lb, that the <CN* > of In is particularly 
reduced by the O-deficiency. It indicates that the O-deficiency should be more easily achieved 
in the presence of In atoms. The result can be understood by the increase in bond strength 
as we go from In-0 to Zn-0 and to Ga-0 [14], consistent with In having the largest atomic 
radius and also with the behavior of the stability of V(O). As shown in Fig. lb and c, both 
AstdCN* and A STD d* around In atoms become all increased by O-deficiency, indicating 
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that as the number of atoms surrounding the large In atoms is reduced, the structural 
distortions around the In atoms becomes enhanced. 

The electronic structures of c-IGZO and a-IGZO are compared in Fig. 2. The localization 
characters of the wave functions at the conduction band minimum (CBM) and valence band 
maximum (VBM) are described by the rod graphs shown in the inset of Fig. 3, which show 
the mapping of wave functions onto atomic sites: we calculated the angular momentum 
decomposition of a wavefunction: \1> = J2i a a i,a(pa, where i denotes atom and a angular 
momentum, and the atomic site mapping q = Yl a a !,«- The site-mappings are parameterized 
by the normalized inverse participation ratio (IPR = (1/N c t )/(V-^ S c ?) 2 ) ( see Fig- 3). 
The IPR increases with the localization and thus can describe the band tail states. The IPR 
of c-IGZO are usually small. 

An important finding is that the fluctuation in the electronic structure induced by the 
disorder is very small in the CBM. The character of the electronic structure around the 
CBM in the disordered a-IGZO is calculated to be amazingly close to that in c-IGZO (see 
Fig. 2 and 3), accordingly the IPR data are similar to that in c-IGZO. From these results, 
we suggest that the localized band tail, which is well-known to be formed around the band 
edges in amorphous materials, can be absent around CBM in the ASMO. Thus the mobility 
of electron carrier in a-IGZO can be comparable to that in c-IGZO. In the high ionic oxides, 
the hybridizations between the cation-s and O-p orbitals are weak. The CBM state is 
characterized almost only by the cation-s-orbitals. Thus the connections between cation-s- 
orbitals at CBM are insensitive to the CNs of cations, unless the dangling bond is formed. 
Therefore the character of CBE is insensitive to the disorder. 

The wavefunction of the CBM is more accumulated around In atoms in both c-IGZO 
and a-IGZO (see the inset of Fig. 3a and 3c), indicating that the channel formed among In 
atoms should be more favorable for the electron transport, thus the higher concentration of 
In should be favored for the high mobility. Since the CBM wavefunction remains sufficiently 
well delocalized in the amorphous structure, an effective mass for the free electron carrier 
can be defined. The calculated effective mass of 0.201 mo in a-IGZO is, thus, very close to 
that in c-IGZO(=0.196 mo in an ab-plane, and 0.213 along the c-axis). 

We note that the effect of the disorder is strong for the valence band. The IPR data 
indicate that the wavefunctions around the valence band (VB) becomes seriously localized, 
as shown in Fig. 3c. There is a region around the VBM where the IPR is very large. These 
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indicate, that the localized band tail states are strongly present at the VBM. It is because 
the VBM orbital consists of strongly directional O-p orbitals. It leads to a lowering of band 
gap. This should be an important feature of the a-IGZO and the ASMO. The band gap 
(Eg) of a-IGZO is calculated to be smaller by 0.44 eV than that of c-IGZO. The localized 
tail states should trap hole carriers (not for electron carriers), thus the low mobility of hole 
carrier should be low. These states can play a role in the illumination stress instability of 
the TFT [3, 7]. 

We now discuss the effect of O-deficiency (one-O-missing in supercell) on the electronic 
states, which is compared to that of c-IGZO. The O-deficiency plays a role of donor. Usually 
the donors creates donor level below the CBM, however surprisingly we could not find any 
localized states around CBM. The atomic mapping and the IPR value of the wavefunctions 
around CBM, is little affected by the O-deficiency. These indicate that a localized orbital is 
not generated by the O-deficiency, i.e. in contrast to the case in c-IGZO. The O-deficiency 
can be described by the change in the substructures around cations rather than by the 
formation of O-vacancy structure, and just as the CBM is insensitive to the disorder char- 
acterized by the admixture of many kinds of substructures, it is not also affected much by 
the change in substructures. Even if a donor-center is not formed, free electron carriers are 
generated due to the chemical imbalance induced by the O-deficiency. 

On the other hand, in c-IGZO, an O-vacancy gives rise to dangling bonds and create the 
localized donor-like trap centers. In the c-IGZO which has a layered structure, four kinds 
of O-vacancy structures can be distinguished. The most stable V(O) is surrounded by three 
In atoms and one Zn atoms. The V(O) preferentially forms around In atoms [15]. This 
trend is consistent with the fact that the O atoms can be easily removed from around In in 
O-deficient a-IGZO. We find that the formation enthalpies of the O-vacancies in c-IGZO are 
much larger that in a-IGZO. The formation enthalpies of these O-vacancies are calculated 
to be 4.14~4.64 eV. It indicates that O-deficiency can be much more easily achieved in 
the a-IGZO than in c-IGZO. Even if the electronic structure calculations indicate that the 
lowest energy V(O) in c-IGZO induces a shallow level, close to CBM, as shown in Fig. 3b, 
the atomic mapping of the wavefunction indicate that the CBM orbital is slightly localized 
at In atoms around V( O), resulting in the increase of the IPR values. This is in contrast to 
the case of O-deficient a-IGZO. This result can provide an explanation for how the carrier 
mobility can be measured to be lower in the c-IGZO rather than in a-IGZO, unless the 
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O-deficiency is too serious [4]. The mobility can increase with the carrier concentration by 
the increasing O-deficiency concentration [4], since the furthermore formation of trap centers 
can be avoided by the O-deficiency in a-IGZO[2]. 

For the valence band, the effect of the O-deficiency is to create the localized states more 
above the VBM. Accordingly the calculated band gap (E g ) was more lowered by 0.59 eV, 
compared to the stoichiometric shown by Fig. 2c. These states can provide an 

explanation for the absorption of blue-green light. As the O-deficiency concentration becomes 
more serious by two-more missing-0 atoms per supercell, the structural fluctuation becomes 
more pronounced and more complex, and sub- structures with weak or dangling bonds are 
formed. Some d* are calculated to be much larger than the fitting line for the a-IGZO with 
the two-more missing atoms, as shown in Fig. la. Both Astd d* and AstdCN* becomes 
distinctly larger around In atoms by the increased O-deficiency (see Fig. 1), indicating 
that the local structures around In atoms become remarkably distorted and thus weak or 
dangling bond tend to be generated around In. In the case of two missing oxygens, we 
identify an asymmetric In04 structure having an isolated dangling bond (see the electron 
density picture in Fig. 3e). It induces a deep donor level within gap, trapping two electrons. 
Furthermore the conduction band minimum is also serious altered, that should affect the 
electron mobility. We repeated several SA processes, but could not avoid the deep level 
formation. The reduction of the electronic energy resulting from the capture of electrons at 
the deep level may stabilize this kind of DB formation. This kind of deep level may lead to 
the Fermi level pinning in the seriously O-deficient a-IGZO, and may give an explanation 
for the subgap states[2], and also reduce the mobility. We note that the subgap states can 
also come from the metastable vacancy structures. 

In conclusion, we investigated the microscopic and electronic structures of a-IGZO as 
a representative of ASMO. We showed that the band tail states, which are well-known to 
be present in amorphous semiconductor, can be absent around the CBM in amorphous 
IGZO, even if the structure is complicated by the admixture of many different kinds of 
substructures. The properties of the CBM are found to be little affected even by the donor- 
like O-deficiency, which is in contrast to the crystalline IGZO, where an O-vacancy forms 
a localized state. These give an explanation to why the mobility of amorphous IGZO can 
be much higher than that of crystalline IGZO. The disorder and O-deficiency are found to 
affect mainly the character of the valence band, resulting in the serious formation of localized 
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states around the VBM. Our results provide a new perspective on the physics of the ASMO. 
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FIG. 1. (Color online) (a) The effective coordination number (CN*) versus the effective bond 
length(d*) in the stoichiometric a-IGZO are shown by solid symbols. The O-deficient a-IGZO with 
one and two missing oxygen atoms are also shown, respectively, by gray and white symbols. The 
fitting lines are obtained by the importance samplings. The large half-filled symbols represent the 
corresponding values in c-IGZO. (b) The average of the coordination numbers (<CN* >) versus 
x, with the standard deviations (AstdCN*) indicated by the error bar, and (c) the standard 
deviations of the calculated d* from the fitting lines (Astd d*) are shown. 
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FIG. 2. The electronic structure of (a) crystalline c-InGaZnO-4, (b) stoichiometric amorphous a- 
InGaZn04 and (c) an O-deficient a-InGaZn04 by one missing oxygen in supercell, where the dotted 
line describes the Fermi level. The Brillouin zone is artificial, given by the rhombohedral supercell, 
and the Z is (0,0, 0.104 A" 1 ). 
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FIG. 3. (Color online) The density of electronic states of (a) the c-IGZO and (b) the O-vacancy in 
c-IGZO, (c) the stoichiometric a-IGZO, (d) an O-deficient a-IGZO with one missing oxygen atom 
(e) an O-deficient a-IGZO with two missing oxygen atoms in supercell are shown. The bar graphs 
in insets indicate the atomic mappings of the wave functions at conduction band minimum (brown 
bars) and valence band maximum (black bars) (see the detail in text) and the localization of the 
wave functions are described by the inverse participation ratio (IPR) calculated from the atomic 
mapping of the wavefunctions (see the text), and the red dotted lines indicate the Fermi level. 
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